arXiv: 1508.02944vl [cond-mat.quant-gas] 12 Aug 2015 


Topological phases of shaken quantum Ising lattices 


Samuel Femandez-Lorenzo,*’ * Juan Jose Garcia-Ripoll,^ and Diego Porras*’ ^ 

^ Department of Physics and Astronomy, University of Sussex, Palmer, Brighton BN19QH, United Kingdom 
^Instituto de Fisica Fundamental IFF-CSIC, Calle Serrano 113b, Madrid E-28039, Spain 

The quantum compass model consists of a two-dimensional square spin lattice where the orientation of the 
spin-spin interactions depends on the spatial direction of the bonds. It has remarkable symmetry properties 
and the ground state shows topological degeneracy. The implementation of the quantum compass model in 
quantum simulation setups like ultracold atoms and trapped ions is far from trivial, since spin interactions in 
those sytems typically are independent of the spatial direction. Ising spin interactions, on the contrary, can be 
induced and controlled in atomic setups with state-of-the art experimental techniques. In this work, we show 
how the quantum compass model on a rectangular lattice can be simulated by the use of the photon-assisted 
tunneling induced by periodic drivings on a quantum Ising spin model. We describe a procedure to adiabatically 
prepare one of the doubly-degenerate ground states of this model by adiabatically ramping down a transverse 
magnetic field, with surprising differences depending on the parity of the lattice size. Exact diagonalizations 
confirm the validity of this approach for small lattices. Specific implementations of this scheme are presented 
with ultracold atoms in optical lattices in the Mott insulator regime, as well as with Rydberg atoms. 


I. INTRODUCTION 

In the pursuit of the quantum computer, the problem of de¬ 
coherence arises as the main obstacle to preserve coherent lin¬ 
ear superpositions as to take advantage of the computational 
power they can provide us with. In principle, quantum er¬ 
ror correction codes offer a solution to achieve a fault-tolerant 
quantum computation [1]. An alternative route consists of us¬ 
ing topologically protected Hilbert spaces [2, 3]. In this con¬ 
text, a two dimensional quantum compass model on a square 
lattice was proposed by Dou^ot et al. [4] as a simple model 
to implement a protected qubit. Generally speaking, ‘compass 
models’ refers to a broad type of lattice Hamiltonians in which 
the couplings between sites depend on the orientation of the 
bonds in the lattice. A thorough review of these Hamiltonians 
and their properties can be found in Ref. [5]. 

The quantum compass model was originally introduced in 
1982 as a toy model to gain insight in the context of Mott 
insulating transition metal compounds, for which one finds 
anisotropy of the exchange for different pairs of ions. The 
name ‘compass model’ arises by analogy with the dipolar cou¬ 
pling in a classical model of magnetic needles arranged in a 
lattice [6]. The 2D-version of this model on a n x m lattice is 
defined by the following spin Hamiltonian (5=1 /2), 

j j 

where CJ;’^ are the usual Pauli matrices and j runs over the 
lattice sites. We shall assume free boundary conditions in the 
following. We can choose ferromagnetic couplings {Jx,Jy > 0) 
without lose of generality, since ferromagnetic and antifer¬ 
romagnetic quantum compass models are related by unitary 
transformations. One observes in Hamiltonian (1) that there 
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are two competing tendencies owing to two types of Ising- 
like interactions; bonds along the y axis induce spin align¬ 
ment along y ((crj) 0), while bonds along the x axis induce 
spin alignment along x ((<Tp 7 ^ 0); the resulting ground state 
is therefore a highly entangled state without an obvious order 
parameter. Many recent numerical studies have examined the 
quantum phase transition of the anisotropic model on a square 
lattice through the isotropic point {J^ = Jy), pointing to the ex¬ 
istence of a first order quantum phase transition [7-9]. It has 
been theoretically shown that this model arises as an effective 
description of the low-energy physics in systems of magneti¬ 
cally frustrated Josephson junction arrays [3, 4], and experi¬ 
ments have shown signatures of the physics of this model in 
few qubit setups [10]. 

The physical implementation of the quantum compass 
model (1) in atomic systems would represent a significant 
breakthrough in research on topologically protected qubits. 
However, experimental tecniques for analogue quantum sim¬ 
ulation [11] typically provide us with effective spin-spin in¬ 
teractions that are independent of the spatial direction of the 
bonds. For instance, Ising interactions, with couplings of the 
form GjOj along every spatial direction, can be readily in¬ 
duced and controlled in systems like trapped ions [12-14], 
ultracold atoms in optical lattices [15-18] and Rydberg atoms 
[19-21]. In this work, we surpass this limitation, by showing 
that the quantum compass interactions can be implemented 
by dressing Ising spin-spin interactions with periodic driving 
fields. The basic idea of our work is the dressing of Ising 
interactions by the photon-assisted tunneling induced by pe¬ 
riodic drivings with a site-dependent phase over a square spin 
lattice . A judicious choice of the site-dependence of the driv¬ 
ing phase leads to the spatial dependence of interactions in the 
quantum compass model, (see Fig. 1). The necessary periodic 
drivings can be implemented with running spin-dependent op¬ 
tical potentials such as those demonstrated in [22]. Our ideas 
have a direct application in the implementation of topological 
models with ultracold atoms and Rydberg atoms with state- 
of-the-art techniques, since basically we only request an ad¬ 
ditional spin-dependent moving lattice to the trapping optical 
lattice potential. Indeed, periodically driven atomic lattices 
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with site dependent phases have brought a lot of attention in 
recent years, as they can be used in the simulation of syn¬ 
thetic gauge fields [23-31]. The dressing of one-dimensional 
quantum Ising systems by periodic drivings with a gradient in 
intensity has been considered recently in [32]. 




Figure 1. General scheme: Our model assumes an initial quan¬ 
tum Ising lattice as a starting point. Periodic drivings with a site- 
dependent phase allow us to dress the original Ising interaction to 
obtain the quantum compass model, in which sites interact through 
spin-components depending on the orientation of the bond. 

This article is organized as follows. In Sec.II we re¬ 
view some general symmetry properties of the quantum com¬ 
pass model that will be necessary to understand the adiabatic 
preparation of the ground state. In Sec.Ill the photon-assisted 
tunneling scheme is applied to engineer the quantum compass 
model using a 2D Ising model together with a convenient pe¬ 
riodic driving. In Sec.IV we discuss an adiabatic preparation 
of the ground state of this model; the corresponding imple¬ 
mentation of this procedure using ultracold atoms or Rydberg 
atoms in optical lattices is presented in Sec.V. Finally, Sec.VI 
summarizes the main results obtained in this work. 


II. SYMMETRY PROPERTIES OE THE COMPASS 
MODEL 


Let us now look over the main symmetry properties of the 
quantum compass model that we shall use throughout this 
work [4]. Generally speaking, a symmetry operation repre¬ 
sented by certain operator P commuting with the Hamiltonian 
H, may lead to the presence of degenerate states. This means 
that P acting on a state 1 1 //) results in another state | <p) with 
the same energy, although this argument fails in the case we 
get the same state, jy/) = |^). Nonetheless, if we find another 
symmetry operation Q such that [R, 2] ^ 0 for any state jy/), 
one can be sure that all states are at least doubly degenerate. 
To see how this works, let us suppose that acting on the state 
I y/) with P and Q we obtain the same state | y/). This implies 
that PQ\^f) = QP\y) = I*/)’ hence [PtQ\ = 0, which is con¬ 
trary to our initial statement. We have therefore proved that 
the resulting states are necessarily different. Leaving aside 
accidental degeneracy, another condition must be imposed if 
we wish to get doubly degenerate states avoiding further de¬ 
generacy. This condition turns out to be [P^,Q\ = [P^Q^] = 0 
[4]. The reason for this is that if one starts with an eigenstate 
of Q{P) and then acting on it with P{Q), the resulting state 


has to be different from the original one as we just proved, 
but acting again on this state with P{Q), on should come back 
to the initial state. Having two sets of non-commuting opera¬ 
tors {P,} and {Qj}, the previous conditions are generalized as 
[Pi, 2;] 7 ^ 0 V/, j and , Q^] = [Pi,Q jQk] = 0 V/, j,k [4] . 

The Hamiltonian (1) has two sets of discrete symmetries 
satisfying the above conditions, namely 
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(2) 

Qh-I\^Ux,jy) 

Jy 

jx = l,2,...,n 

(3) 


i.e, each Pj^ is the row product of CJ' in that row while 
is the column product of C7' in the column. Physically, 
Qj^ corresponds to a rotation by an angle K around the x 
axis of all the spins of the column labeled by jx, while Pjy 
corresponds to a rotation by an angle n around the y axis 
of all the spins of the row labeled by jy. In particular, they 
satisfy {Pjy,Qj^} = 0 Wjxjy and [PiPj,Qk] = [Pi,QjQk] = 0 
yi,j,k. Thus, aside from accidental degeneracies, one expects 
every state to be doubly and only doubly degenerate. We 
shall then assume that the ground state of the quantum 
compass model is effectively a two level system satisfying 
the conditions of a protected qubit. Local noise acting on 
a single lattice site may not commute with the symmetry 
operators Pjy and Qj^ corresponding to the row and column 
of such site, but the remaining symmetries ensures the system 
to remain doubly degenerate. This is the case unless at least 
min(OT,n) of these perturbations act simultaneously over 
the whole lattice. The set of integrals of motion Pj^ (or 
alternatively Qj^) can be used to distinguish between the two 
degenerate states of the ground state since they have different 
quantum numbers, namely either (pi,..., p„,) = (1,..., 1) 
or (pi,...,pm) = This result was proved in 

Ref [33] for a square lattice, and it can be straightforwardly 
generalized for a rectangular lattice. Heuristically, one may 
expect this result by exploring the trivial case for which 
Jx = 0. In such a case the model is simplified to a set of Ising 
columns with ferromagnetic coupling. For simplicity, taking 
the square lattice Ixl, the ground state consists of 2' states 
defined by = ... = mj^i = ±1 with jy = 1,...,/, where 
nij is the eigenvalue of dj. As all the rows are identical, all 
the are equal as well for all jy, thus we conclude the 
value of Pjy is either -|-1 or —1 for every row. 


III. PHOTON-ASSISTED TUNNELING 

In this section we shall use the photon-assisted tunneling 
toolbox [24] to implement the quantum compass model. To 
understand how to achieve this goal, let us first have a look at 
the quantum Ising model. 

Ho = 

j j 


(4) 
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Expressing the spin operators in terms of the raising and low¬ 
ering operators, i.e. = (<7+ -t-cr^) and a-' = —i((7+ — <7^), 
Eq. (4) turns out to be, 

Ho = -J. E ^j+e) 

J 

-JyJl {<^1 <yj+ey + ( 5 ) 

J 

In contrast, rewriting the quantum compass model (1) in a 
similar way, we obtain a slightly different expression H.c. 

= - E cTf^e.+Jt<yt <yi+e,) 

j 

-L{-Jr^j^!+ey+Jr<yj<yj+ey)+ii-^-^ ( 6 ) 

j 

where we have defined /++ = 7+^“ = Jx,y- By the use of the 
photon assisted tunneling in the original Hamiltonian (4), we 
aim for finding a set of effective coupling constants such that 
Eq.(5) equals Eq. (6), which in turn means 

(/7+)eff=(4+^)eff (7) 

(77^)eff = —{Jy )eff- 

Now that our goal is clear, let us use the ingredients of the 
photon-assisted tunneling toolbox to find how to achieve the 
conditions (7). In doing so, the Hamiltonian of the system can 
be written as, 

H = Ho+HdiT), (8) 

where we define, 

Hd{T:) = Y,^^j+lli^dCos{cooT + ^j)a'j. (9) 

j j 

The second term in (9) represents a periodic energy driving 
of the qubit, while Clj is typically chosen such that the first 
term in Eq.(9) represents a gradient of the individual frequen¬ 
cies [23], i.e., Q.j = Ho -f • j, although it will be taken as 
a constant for the purpose of this work, Clj = Q.. These two 
elements are all that we need to take advantage of the photon- 
assisted tunneling toolbox [24]. Notice that there is a freedom 
in choosing the spatial dependence in (j)j, and we shall assume 
a linear dependence in both x and y axis, 

(l)j ^A<j)xjx+A(j)yjy, (10) 

where A(l)x^y are given constants, and jxj are positions in the 
lattice. 

Let us now express the Hamiltonian Hq in the interac¬ 
tion picture with respect to //j, namely Ho{t) = U[t)^HoU (f) 
where U{t) = In this picture, the raising and 

lowering operators evolve like. 


Notice that the last term in Eq.(ll) can be gauged away us¬ 
ing the unitary transformation (7^ —jqgjjj-g Qjjg 
obtains a Hamiltonian having the same structure as Eq.(6), in 
which we replace the bare couplings by their corresponding 
time-dependent dressed couplings 

(4V)eff = W). (13) 

To help in the analytical treatment, we make use of the Jacobi- 
Anger expansion, 

(14) 

—oo 

where are Bessel functions of the first kind, yielding 

s,s' 

(15) 

{JtyU = Jx,yY^Js{n)JAn)e'^^’^^^+^^^e-^ 

s,s' 

(16) 

Assuming Jx^y H and tuning the driving frequency to (Od = 
2Q., the rotating-wave approximation (RWA) allows us to ne¬ 
glect the fast-oscillating terms and to keep those terms fulfill¬ 
ing the resonance condition = i -|- 1 in Eq.(15) and s' = s 
in Eq.(16). Having done this, we obtain the following dressed 
couplings, 

<h^^J+exy\ 

{J++U = Jx,y^++{T]A^x,y)e V ^ ) (17) 

(77- )eff = Jx,y^^,; (I?, A^x,y). ( 18 ) 

where we define a set of complex amplitudes, 

,^7,+ (77,A<^,,,) (19) 

S 

^x,y{nA^x,y) = E</7’?)e^'^- (20) 

A 

Einally, rewriting the term ^g j^ay use 

the following gauge transformation: fft i —y cjtresult¬ 
ing in the effective Hamiltonian we were looking for, 

Heff =-JxE V.0 

J 

-7,E +H.C. (21) 

J 

Eollowing the idea given in Eq.(7), we notice that the quantum 
compass model can be implemented if we manage to find a set 
of parameters, (rj ,A^x,A^y), such that 

^x^+(il,A(l>x)=^x^-(V,A(^x,) (22) 

^++iri,A(j)y) = -^+-{ri,A(^y). 
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Is there any solution to the system of equations (22)? We 
will show numerically that there are actually an infinite num¬ 
ber of solutions. First, note that Eq. (22) can be expressed 
as, 

I i'n , ^(l>x,y) I = I I *^23) 

arg(,^/+(T7,A0;c)) = arg(,^/^(77,A0^)) -f 2^;r (24) 

arg(,^++(T7,A0v)) = arg(,^+^(T7,A0y)) -f ^;r. (25) 

To solve this system of equations, let us define the func¬ 
tions f(j], A(j)^^y) = I ^+y+ (j] , A0^) I / I ^+y- (77 , A0^,y) I - 1 and 
g{nA^x,y) = srg{^+y+{ri,A(^y))-arg{^+y-{ri,A(l)y)). Triv¬ 
ially, we see that Eq.(23) is equivalent to find the zeros of 
the function /, while the solutions to Eq.(24) and Eq.(25) are 
equivalent to the function g taking the values kn and 2kn re¬ 
spectively. The jellyfish-like pattern showed in figure 2(a) rep¬ 
resents the function /, for which we have limited the range of 
the function to the interval [—0.1,0.1] so that the graph clearly 
shows the region in which the function becomes zero. On 
the other hand, the function g is shown in Fig.2(b), limited 
to the range [0,2;r]. In looking at the graphs, it is straight¬ 
forward to confirm that there are an infinite number of so¬ 
lutions. The isotropic point, (Jx = Jy) is reached when the 
solutions are taken symmetrically with respect to the sym¬ 
metry axis of the jellyfish-like pattern; an example could be 
{t] = l,A^x = 4.6539, A(j)j. = 1.6293}, which gives the value 
,^++ = ^+- = 0.5368. 

IV. ADIABATIC PASSAGE 

Now that we have shown how the quantum compass model 
can be implemented using periodic drivings, we aim for 
preparing the ground state of this model. One way to do 
so is by finding an adequate adiabatic passage starting from 
the ground state of certain initial Hamiltonian for which the 
ground state is known and can be prepared, and then slowly 
changing this Hamiltonian until the quantum compass model 
is eventually reached. One possible option using a transverse 
magnetic field will be discussed in this section. 

Let us then consider an additional transverse magnetic field 
along z in our original Hamiltonian (8) as follows, 

H^Ho+Hj{T) + '£^a], (26) 

j 

where 5 is a time-dependent parameter that measures the 
strength of the field. Notice that we had already introduced 
a transverse field along z in Eq.(8), so the overall transverse 
field depends on the sum of both terms, e/2 — {Q. + 5)/2. One 
may also express Hq -f 5(7} in the interaction picture with 
respect Hd', it is then straightforward to check that, following 
the same procedure used in Sec.Ill and setting the same reso¬ 
nance condition (Od — 2fi, we arrive at the following effective 
Hamiltonian, 

^eff(0=^^eff + i:^C7j, (27) 

j 


f(T1,A(|)) 



U 

(b) 

Figure 2. In this figure we plot the functions that govern the spin- 
spin interactions as a function of the periodic driving parameters, 
(a) Function /(p, A(/).i,j,) = |,^++(p,A(|)^)| / (p,A(/)x, 3 ,)| - 1. 

The values f{ri,A(j>xj) = 0 represents the solutions to Eq.(23). 
The range of the function is limited to the interval [—0.1,0.!] to 
show clearly the points at which the funcion is zero, (b) Func¬ 
tion g(ri,A(j>xy) = arg(,^++ (p,A((),,)) - aig(,^+f{ri,A<j)y)), limited 
to the range [0,2;r]. The regions with g{ri,A(j)xy) = 0 correspond to 
the solutions to Eq.(24) and Eq.(25). 


where was given in Eq.(21). Therefore, ensuring the 
magnetic field is strong enough so that 5 >> H, p >> Jx,Jy, 
one can safely assume that the magnetic field is the domi¬ 
nant term in (27), being the ground state |0) = Hy I t)y in 
terms of eigenstates of CTj. Finally, a feasible adiabatic pas¬ 
sage to reach the compass model consists of decreasing e from 
£ >> Q, 7] slowly enough until reaching the compass condi¬ 
tion z — £l — (i>dl2. (5 = 0). One expects this system to un¬ 
dergo a quantum phase transition as the magnetic order of the 
transverse magnetic field and the quantum compass model are 
different. 

The adiabatic approximation describes, upon certain condi¬ 
tions, how slowly we need to vary 5(f) to ensure the system 
remains in the ground state through the evolution. Accord¬ 
ing to this approximation, the instantaneous eigenstates of the 
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time-dependent Hamiltonian H{t) at a given time evolve con¬ 
tinuously to the corresponding eigenstates at later times, pro¬ 
vided that the eigenenergies do not cross and the evolution 
is slow enough. The intrinsic time scale used to determine 
what slow and fast mean is usually provided by the gaps in 
the spectrum. This also provides a general validity condition 
for adiabatic behaviour that corresponds to the probability that 
the final state of the system is different from the initial state 
[34], 


max 

0 <t<T 


{k\Hit)\n) 

8nk 


« min 

o<f<r 


(28) 


where T is the total evolution time, and g„i^ the energy gap 
between level n and k. Note that the ground state of the trans¬ 
verse field Hamiltonian is unique, while the compass ground 
state is doubly-degenerate, thus we expect crossing of levels 
at the end of the adiabatic preparation; moreover, the natural 
question arises about what state is eventually reached. 

A simple possible way to overcome the issues above is by 
hnding an integral of motion conserved through the evolution 
that we may use to tell the twofold ground state apart. Un¬ 
der this assumption, the adiabatic approximation would still 
be valid if there is no crossing of levels with the same con¬ 
served quantum number. As it was shown above, the set of 
integrals of motion Pjy{Qj^) can be used to distinguish be¬ 
tween the two degenerate states, since they have different 
quantum numbers, namely either (pi,... ,pm) = (1, • • •, 1) or 
(pi,..., Pm) = (— 1 , ■ ■ • 5 — 1 ) (similarly for the quantum num¬ 
bers {qi,... ; however these are not good quantum num¬ 

bers when the magnetic field term is included in the Hamil¬ 
tonian. The operator Z = is fortunately an integral of 

motion of the whole Hamiltonian (27) as it is straightforward 
to check, and so it can be used in principle to determine the 
actual ground state reached at the end of the adiabatic passage. 
Surprisingly, the parity operator Z can distinguish between the 
twofold ground state only in those cases such that we have an 
odd-odd or odd-even lattice size. In Appendix I, it is proved 
that for both odd-odd and odd-even cases the eigenvalues of Z 
are either z = -|-l orz = — 1 , while in the even-even case the 
eigenvalue is always z = +1. Hence we expect two different 
kinds of behaviour in the adiabatic passage depending on the 
parity of the lattice size. In the hrst case, the two ground states 
of the compass model are not connected, thereby the only rele¬ 
vant gap to be considered is between the ground and hrst even 
excited state. Recall that the ground state of the transverse 
held Hamiltonian has eigenvalue z = +1, so we expect to pre¬ 
pare the ground state of the compass model corresponding to 
this eigenvalue; however we would be unable to prepare the 
ground state corresponding to the eigenvalue z = — 1 using this 
approach. In the second case, with an even-even lattice size, 
any superposition of the ground states are eigenstates of the 
operator Z with eigenvalue z = -fl, and for that reason, in this 
case we could only assume that the hnal state is a superposi¬ 
tion state given by the specihc adiabatic evolution performed 
on the system. The adibaticity of this case would be given by 
the gap between the hrst and and second even excited states. 

Exact diagonalizations for small systems (up to a 5x4 lat¬ 
tice) were performed to observe the evolution of the relevant 
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Figure 3. Quantum phase transition from transverse held to compass 
model for (a) even-even lattice size and (b) odd-even/odd-odd lattice 
size. We show the gaps together with the absolute value of the nor¬ 
malized magnetization along z of the model (29) for different values 
of the lattice size. The blue and green lines corresponds, respectively, 
to the gap between the hrst and second even excited states with re¬ 
spect to the ground state. The gaps in Figure (b) are represented in 
a semi-log plot. The value A = 1 corresponds to the initial trans¬ 
verse held Hamiltonian and A = 0 represents the isotropic quantum 
compass model. 
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gaps involved in the quantum phase transition so as to exam¬ 
ine the validity of the adiabatic approximation. Additionally, 
the magnetization along z was found to be an order parame¬ 
ter of the quantum phase transition, where the ground state of 
the transverse magnetic field is magnetically ordered, ^ 0 , 
and the compass ground state is disordered, = 0 . In order 
to simplify the analysis, the following parametrization was in¬ 
troduced, 

// = (l-A)i/c + A^m (29) 

f j 

Hm — E ’ 

j 

therefore A = 1 corresponds to the initial transverse field and 
A = 0 represents the isotropic quantum compass model. 

Keeping in mind the results for the integral of motion Z, 
the relevant gap in an odd-odd and odd-even lattice is the one 
given by the first even excited state of the system, A. In con¬ 
trast, this gap goes to zero in an even-even lattice, so the rele¬ 
vant gap in this case. A', is given by the energy difference with 
the second even excited state. Note that those levels with odd 
parity (z = — 1) do not play any role in any case as the Hamil¬ 
tonian does not connect states with different parities. Fig. 3 
shows the evolution of these gaps for different sizes together 
with the magnetization \M^\ as the parameter A changes from 
A = 1 to A = 0. For the odd-odd and odd-even cases we con¬ 
firm that there exists a finite gap along the adiabatic passage 
between the ground state and the first even excited state that 
narrows as the lattice size becomes larger. As expected, this 
gap goes to zero in the even-even case when the two degener¬ 
ate even ground states collide in the compass model as shown 
in Fig. 3(b). Therefore, the relevant gap when we approach 
A = 0 is given by the energy difference between the first and 
second even excited states A'. In concluasion, these results 
show that the even ground state of the compass model can 
be prepared using this adiabatic passage in the odd-odd and 
odd-even cases , whereas for the even-even case we expect 
the system to be in a superposition of the ground states given 
by the specific evolution of 5(f). 

V. PHYSICAL IMPLEMENTATIONS 

In this section we show how our ideas can be implemented 
using specific atomic experimental setups. One possible im¬ 
plementation of the quantum compass model using Josephson 
junctions arrays was proposed by D 0 u 50 t.et.al [4]. A close re¬ 
lated model was actually implemented in a proof-of-principle 
experiment using superconducting nanocircuits. 

Atomic systems present many advantages for quantum state 
preparation and measurement. Thus, a scalable and efficient 
implementation of the quantum compass model in atomic ex¬ 
perimental setups would be very useful. Furthermore, im¬ 
plementing a controllable longitudinal magnetic field like in 
Eq. (30) may allow experimentalist to adiabatically create the 
topologically degenerate ground state. Periodic drivings like 


those required for our proposal can be implemented in atomic 
systems by means of lasers, with a site dependent phase that 
corresponds to the laser optical phase. Ions trapped in two- 
dimensional arrays of microtraps or Coulomb crystals could 
be considered here, because optical forces can be used to in¬ 
duce Ising interactions [12, 14]. However, the dipolar decay 
of trapped ion spin-spin interactions would lead to long-range 
quantum compass models, with properties that may depart 
from the original Hamiltonian (1). In the following we fo¬ 
cus on proposals for atomic setups that may provide us with 
short-range Ising interactions, namely, neutral bosonic ultra¬ 
cold atoms and Rydberg atoms in optical lattices. 

A. Ultracold bosons in optical lattices 

We need first to understand how the following effective 
quantum Ising Hamiltonian can be implemented with ultra¬ 
cold bosons. For this we rely on a quantum simulation pro¬ 
posal which relies on the internal state of atoms that are frozen 
in a Mott insulator state in an optical lattice [15, 16]. Un¬ 
der appropriate circumstances, hopping can only be a virtual 
process that enables superexchange interactions, as explained 
in [16] through perturbative calculations. These interactions 
have been demonstrated experimentally, both in superlattices 
[35, 36], as well as in longer tubes with a few quasiparticle 
excitations [37, 38]. 

We consider an optical lattice in the Mott-insulator regime 
with one atom per site (unity filling) and each atom having 
two accessible internal states, playing the role of a pseudo¬ 
spin S = 1/2. The atoms may be formally identified with two 
types of bosons,' one may denote the bosonic 

operators a,- and 5,- as the destruction operators of each inter¬ 
nal state at the site i. Such a system is well described with 
the Bose-Hubbard Hamiltonian when the energies involved 
are small enough so that the second Bloch band never gets 
populated, 

^=E^f7 + E^f'' (30) 

(‘J> j 

= —Ja{a\aj+h.c.)+Jb{b\bj + h.c.) 

^ \Uaaa]a^ja jOy + ^ Uhhb'^jb'^jbjbj + Uaba^jb]bjaj, 

The Hamiltonians and represent, respectively, the 
probability of atoms hopping to neighbouring sites and their 
effective on-site interaction. Assuming we are in the Mott- 
insulator regime, 7 <C U, the hopping Hamiltonian can be 
considered as a small perturbation with respect to all other 
terms. We project our problem into the subspace of sin¬ 
gle atomic occupation and use quasi-degenerate second-order 
perturbation theory to obtain an effective spin Hamiltonian. 
The corresponding operators of the effective spin system are 
< 7 / = a\bj, aj = ajb\ and < 7 ) = a\aj — b\b;. The effective 

J J J ■'7 J J J •' 

spin Hamiltonian reads 

= E (A"(7/C7j -f A-L(<7/(7/ -f -f (31) 

{‘J) • 
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with constants given by [16] 


1 =^ = 


il I j2 

■'a 

2Uab 


b _ ^ _j ± _ 

Uaa Ubb 


-2 


JaJb 

~u7b' 


Uaa Ubb 


(32) 


For details on the range of validity and the derivation of those 
equations we refer the reader to Ref. [ 16] . We notice that if the 
term = 0, then we can use the remaining Ising interaction 
as a starting point to derive the quantum compass model. To 
be in a regime such = 0, one could suppress the tunneling 
of h-atoms by choosing a spin-dependent lattice with Jb^Ja- 
If we neglect Jb we get 

= + (33) 

('■7> '■ 


with J = Ja {I/Uaa “ l/(2Uab)), and h~ = -Jl/Uaa- 

To implement periodic drivings, we consider additional 
lasers inducing Raman transitions between levels f and ],. 
Here, the spatial dependence of the phase will appear natu¬ 
rally, since the optical phase of the lasers vary linearly from 
site to site. We need to implement the driving and longitudinal 
field in Eq. (9). Note that, since we obtained an Ising interac¬ 
tion in the z-basis, the longitudinal helds required in Eq. (9) 
must be expressed in terms of <T* operators. A constant field is 
implemented by a two-photon Raman transition or microwave 
field inducing transitions between the two atomic levels, 

= (34) 


The periodic driving fields are then implemented by a 
running-wave potential induced by pairs of lasers with effec¬ 
tive wavevector Ak, and relative detuning within each pair Od, 

Mas(f) = (35) 


Note that time-dependent optical lattice potentials have been 
implemented for spin-dependent transport of ultracold bosons 
in optical lattices (see for example [22]). We notice that the 
laser optical phase is translated into a site-dependent optical 
phase with a linear dependence on the site position [23, 24], 

0,- = ufo(Ak/), (36) 


where we have used that r,- = d^i, with do, the distance between 
sites in the lattice. After implementing a rotation of the spin 
basis (a^ = CJ^, a' = —CJ^ ) in Hamiltonians (33, 34, 35), we 
obtain the driving term in Eq. (9), with qWd = fiias, and the 
phase gradient A(j)xj = {Ak)xjdo. Our final Hamiltonian is 

H = //s +7/mw +77las(0 = 

- E +E^"^f + ? E + 

{i.j) i ^ i 

^5Icos(Akr/-a)df)a?. (37) 


We get our target Hamiltonian plus an additional magnetic 
field term, )■ Under the resonance 

condition Wd = 20 raising and lowering operators rotate with 
frequency Q. (see Eq. (11)). Thus, the term can be ne¬ 
glected in a rotating wave approximation in the limit fi ^ 
J,hA 


B. Rydberg atoms in optical lattices 


Rydberg atoms offer us another physical setup with Ising 
interactions that can be controlled by periodic driving helds. 
We consider a square lattice with one single Rydberg atom 
per site. An effective spin is formed with the states |— )j = 
|g)j, and |-|-)j = |r)j, corresponding to the ground and excited 
Rydberg state, respectively. The Hamiltonian describing this 
system is given by [19] 


Uy Of+lC75+l 

2 ^ 2 2 
<ij> 



(38) 


We have assumed fast decaying interactions between Rydberg 
states, such that the effective spin-spin interaction runs over 
hrst neighbours only. To obtain the compass model we need a 
periodic driving in the jr-basis, since Ising interactions appear 
in the z-basis. Eurthermore, we need to counteract the local 
longitudinal held given by =A + U/2. Eorthis, we consider 
two counter-propagating driving helds inducing a Raman tran¬ 
sition with wavectors k and —k, and detunings (0o + (0d and 
(Oo — (Od, respectively. 

HiUt) = 

I 

^ ^ -PH.c.) . (39) 

We choose C0o = A + U/2to counteract the effect of the longi¬ 
tudinal held. An additional coupling with vanishing effective 
wavevector is required to implement the term proportional to 
fi in Eq. (9). Note that transitions between atomic ground and 
Rydberg states usually required two-photon Raman processes, 
such that the effective wavevector can vanish with a suitable 
orientation of the individual laser beams. We choose a term of 
the form, 

(40) 


If we express the sum of contributions (38,39,40) in a frame 
rotating with the frequency 0 )), we get 


Hrj + Hi 

U 
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Hq, 




<i,i> 


After a rotation of the spin basis we obtain our targeted driven 
Ising model, with J = —U/A, and the phase gradient given by 
the laser optical phases, 0,- = £/o(k i). 
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VI. CONCLUSSIONS 

We have shown that Ising interactions in a square lattice 
can be dressed by a periodic driving field and transformed 
into a quantum compass model. The key idea is to use site- 
dependent driving phases such that the dressed spin-spin inter¬ 
actions depend on the orientation of the bonds connecting lat¬ 
tice sites. We have also show that the ground state of the quan¬ 
tum compass model can be reached by adiabatically ramping 
down a field in the spin z-direction. By using symmetry ar¬ 
guments, we have found conditions under which one of the 
degenerate ground states can be reached, depending on the 
initial quantum state, and the number of sites in the lattice. 

We have discussed two possible implementations with ul¬ 
tracold bosonic atoms and Rydberg atoms in optical lattices. 
However, our ideas can be used in other experimental setups, 
for example in two-dimensional arrays of trapped ions. The 
latter system requires further investigation, since Ising inter¬ 
actions in trapped ion setups are long-ranged [ 12 ], something 
that would lead to the implementation of long-range quantum 
compass models. Any experimental setup with spin interac¬ 
tions where couplings can be dressed by periodic fields is also 
amenable for the implementations of our ideas, like for exam¬ 
ple, arrays of superconducting qubits interacting with classical 
fields in transmission lines [39, 40]. 
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Appendix A: 

Eigenvalues and eigenstates of the operator Z 

In this section we find the eigenvalues and eigenstates of 
the operator Z = H j <7? corresponding to the Hilbert subspace 
spanned by the doubly-degenerate ground state of the quan¬ 
tum compass model. First, recall that the set of integrals of 
motion {Pjy} defined in Eq.(2) can be used to characterize the 
twofold ground state since they have different quantum num¬ 
bers, namely either (pi,... = (1,..., 1,..., 1) or 

-1). If we define the 
states corresponding to the previous eigenvalues as |±)p, this 
result can be written as = ±|±)p V jy. Alternatively, 

the set {Qj^} can also be used for this purpose, giving the 
quantum numbers qj^. The corresponding eigenstates satisfy 
QjJi±)q = ±|±)^ yjx, where the states |±)^ are in general 
different from |±)p. This result was proved in Ref.[33] for a 
square lattice I x I, and it can be straightforwardly generalized 
for a n X m lattice. 


Let us define the operators X = O’! and Y = ]”[;• We 

can express X and Y in terms of the sets of operators {Pjy} and 
{QjA as X = UlQjx and Y = 11% Pjy Thus the states |±)^ 
and |±)p are also eigenstates of X and Y respectively, leading 
to the eigenequations 

A|±), = (±l)”|±), (Al) 

T|±)p = (±)'”|±)p. (A2) 

Notice, though, that the action of X on the basis |±)p is 
unknown, and similarly for Y on the basis |±)^ . Bearing in 
mind the relation (a = 1,2,3) for 

the Pauli matrices, we find a useful equation to relate the pre¬ 
vious operators, XY = This allows us to calculate 

the commutators and anti-commutators for these operators, 
namely 

{X,Y}^2Re{i^^'')Z (A3) 

[X,Y]=2Im{f‘^'')Z. (A4) 

We will show in the following how the above equations en¬ 
able us to determine the eigenstates and eigenvalues of the op¬ 
erator Z in three different cases, depending on the parity of the 
lattice size, i.e. the parity of n and m. Before doing so, let us 
prove a result that we shall require later. We would like to re¬ 
late the basis |±)p and |±)^, and this can be done by the fact 
that {Pjy,Qj,} = 0 yjxjy Effectively, as P;J±)p = ±|±)p 
yjy, this implies Pjy{Qj^\±)p) = =f( 2;J±)/>)- Hence we in¬ 
fer that the state QjJ^±)p has to be proportional to one of the 
eigenstates of Pjy. The only possible solution in this particu¬ 
lar case is Qj^\±)p |T)p, since the choice Qj^\±)p ox |±)p 

leads to contradiction Pjy\±)p = =F|±)p. Using an analogous 
argument for Pjy, one can show that Pjy\±)q ox |=p)^. 

1. Case 1: even n-even m 

In this case the previous relations take the following form. 


{X,Y)=2Z (A5) 

[A,Z] = [y,Z]=0 (A 6 ) 

^|±), = |±), (A7) 

T|±)p = |±)p. (A 8 ) 


Eq.(A 6 ) and Eq.(A 8 ) imply Y (Z|±)p) = Z|±)p. Hence the 
states Z|±)p have to be proportional to the eigenstates of Y, ei¬ 
ther Z|±)p |±)p orZ|±)p |t)p- The proper option can be 

inferred from the anti-commutator relation (A5), which states 
that Z = XY. Effectively, introducing the state |±)p to both 
sides of Z = XY, we find that Z|±)p = A'|±)p. For even n, 
the operator X can be expressed as an even product of Q'^^s 
since X = Ga’ sach of these terms acting in a way that 
Gy^|±)p |=f)p as proved before. Therefore, the only pos¬ 
sibility is that Z|±)p |±)p. The eigenvalues z± can also 
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be calculated using Eq.(A5), since we have p(±|Z|±)p =p 
(±|ZF|±)p, which leads to p(±|±)pZ± =q (±|±)p. Finally, 
we conclude that the states |dz)p are degenerate eigenstates of 
the operator Z with eigenvalues z± = 1. Similarly, the same 
result holds for the the states |±)p, i.e. Z|±)p = |±)p . 


3. Case 3: odd n-odd m 

Finally, in this case our initial relations lead to a familiar set 
of equations. 


2. Case 2: even n-odd m 

The basic equations for this case can be written as follows. 


{Z,T} = 2Z(-1)("/2'''") (A9) 

[X,Z] = [Y,Z]=0 (AlO) 

X\±)q=\±)q (All) 

T|±)p = ±|±)p. (A12) 

Equivalently to the previous case, as we have an even 

product of Q’j^s we infer that Z|±)p |±)p- Now 


using Eq.(A9) the eigenvalues z± can be computed as 
p(±|±)p(-l)("/^^'")z± = (±l)p(±|±)p. Therefore, the oper¬ 
ator Z satisfies the relation Z|±)p = ±(—l)("/^^"’^|±)p. This 
case is trivially equivalent to the odd-even case by a proper ro¬ 
tation, and in such a case the operator Z fulfills an analogous 
relation for the states |±)p. 


{X,Y} = {Y,Z} = {Z,X}=0 

(A13) 

[X,T] =2i“^"Z 

(A14) 

[T,Z] = 

(A15) 

[Z,XJ = 2i“^"T 

(Alb) 

X|±)^ = ±|±)^r 

(A17) 

T|±)p = ±|±)p. 

(A18) 


Notice that these are essentially the same commutation and 
anti-commutation relations as the ones for the Pauli matrices, 
so we can use the well-known results of this representation to 
assure that the eigenvalues Z in this subspace are z± = ±1. 
One can also take advantange of the results for the Pauli ma¬ 
trices to express the eigenvectors of Z, Izb)^ in terms of the 
eigenvectors of X or Y. In particular, in term of the basis |±)p 
we have |±)j = l/^/2{\+)q + \-)q). 
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